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PREFACE 

In  the  preparation  of  this  brief  outHne  of  accepted  methods  of 
dealing  with  Exterior  Ballistics,  the  works  of  Euler,  Otto,  Siacci, 
Parodi,  Ingalls,  Ronca,  and  Charbonnier  have  been  consulted. 
On  account  of  the  finished  classical  niethods  current  at  the  pres- 
ent time  little  remained  for  the  writer  to  do  beyond  selecting 
the  best  for  use  in  this  school.  For  the  appendices  the  writer 
assumes  entire  responsibility.  In  Appendix  I  the  Report  of  the 
Chief  of  the  Weather  Bureau,  for  1900,  and  the  Annuaire  du 
Bureau  des  Longitudes  (1900)  are  used  as  authorities  on  the  tem- 
perature gradient. 

Appendix  II  was  indebted  for  its  basic  hypothesis  to  the  rule 
used  in  France  of  having  the  characteristic  index  indentical  with 
the  sine  of  the  semi-ogival  angle.  The  courtesy  of  Colonel  Rogers 
Birnie,  Ordnance  Department,  U.  S.  Army,  and  of  Commander 
R.  H.  Jackson,  U.S.  Navy,  has  rendered  possible  the  experimental 
verification  of  the  formulas  deduced. 

Appendix  III  is  influenced  partly  by  the  writer's  own  experi- 
ence and  partly  by  the  work  of  Major  Chapel.  The  values  of  the 
function  of  retardation  as  determined  by  the  Gavre  Commission 
were,  of  course,  an  indispensable  element  of  the  discussion. 


ALSTON  Hamilton. 


Fort  Monroe,  Va.,  September  9,  1908. 
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ANALYTICAL 
EXTERIOR  BALLISTICS. 

Satisfactory  analytical  treatment  of  the  flight  of  projectiles 
requires  the  separate  discussion  of  the  principal  problem,  as 
distinguished  from  the  various  secondary  problems  necessary 
for  a  complete  practical  solution. 

The  principal  problem  is  solved  by  using  certain  hypotheses 
which  are  broadly  true  and  which  permit  of  a  discussion  of  the 
problem  without  the  restraint  that  would  ensue  if  all  the  less 
significant  features  of  the  problem  were  introduced. 

The  secondary  problems  are  those  in  which  are  determined 
the  effect  of  the  various  circumstances  omitted  in  the  discussion 
of  the  principal  problem,  in  order  that  proper  compensation 
may  be  made. 

In  Part  I  these  effects  have  been  given  and  also  the  method 
of  compensation  in  each  case.  The  discussion  of  the  secon- 
dary problems  has  been  in  some  instances  omitted  from  Part  I, 
to  be  included  in  Part  II.  * 

In  Part  I.  the  principal  problem  was  hot  solved,  but  the 
tables  resulting  from  its  solution  in  the  cases  of  practical  im- 
portance were  assumed  and  used  in  the  solution  of  ordinary 
problems. 

In  Part  11.  the  principal  problem  is  first  discussed  in  its 
most  general  form  and  the  method  of  solution  outlined. 

Particular  hypotheses  are  then  made  permitting  of  solution 
of  the  problem  for  direct,  curved,  and  high  angle  fire,  and  the 
calculation  of  tables  for  use  in  the  solution  of  ordinary  problems. 

The  Principal  Problem. 
The  principal   problem  in  exterior  ballistics  considers  the 
path  and  circumstances  of  motion  of  a  projectile  moving  with 
a  velocity  v,  at  an  inclination  6  with  the  horizontal,  the  retar- 
dation being  expressed  by     ^~  ,  in  which  F(v)  is  an  increasing 

function  of  v,  and  C  (see  Chapter  II,  Part  I)  an  increasing 
function  of  the  ability  of  the  projectile  to  overcome  the  resis- 
tance of  the  air.  Gravity  is  supposed  constant  in  intensity  and 
acting  vertically  downward  and  the  earth  plane  and  motionless. 
The  air  is  supposed  to  be  of  uniform  density. 
*  See  Appendix  I,  Part  II. 


Under  these  conditions 

dv__  F(v)  "  ... 

dt  C  ^  ^ 

Taking  as  an  origin  some  particular  point  on  the  path  of 
the  projectile  for  which  values  of  v  and  0  will,  for  purposes  of 
distinction,  be  represented  by  V  and  ^p  respectively;  and  repre- 
senting distances  measured  horizontally  in  the  direction  of 
motion  by  x  and  those  measured  in  an  upward  vertical  direction 
by  y,  we  have,  at  any  time,  t  seconds  after  the  projectile  has 
left  the  origin,  for  the  horizontal  and  vertical  components, 
respectively,  of  the  acceleration 

d-'x_d(vcos^)  _       F(v) 


dt-'  dt  C 

dT^d(vsin^)  ^  __  F(v) 
df  dt  C 


cos  0  (2) 

sin  ^  -  g  (3) 


or 


(6) 


COS  ^  ji.    —  V  sm  0  ^,    = ^'  cos  0  (2  ) 

sm  0  ^^  +  V  cos  ^^  ^^  —  —      Q     sm  0  -  g  (3  ) 

From  (2')  and  (3')  we  find  by  combination 

dv  ,   F(v)_     ,      .dO_  ^AO         g  ,.^ 

H \,     —  V  tan  0      —  —V  cot  0^^ .-^-„  (4) 

dt  C  dt  dt       sm  ^  ^  ^ 

whence 

-^  =-  —  g  cos  ^^  (5) 

From  (2  and  (5) ,  on  eliminating  dt 

d  (vcos6')  ^  vF(v) 
d^  gC 

The  accelerations  involved  in  equations  (1),  (2),  (3)  and  (5) 
are  illustrated  in  Figure  1. 

In  the  figure 
PA  =  -  g;  PB  =  —  -^"^^^  ;  PC  the  resultant  of  PA  and  PB. 

PD  =  PB  cos  0  =  ^-^-^-  =^  =  -  ^  cos  0 
dt  dt  C 


From  equation  (5) 


dt 


S 


sec  d  dO 


and  since  dx 


V  cos  (^  dt;  dy  =  V  sin  0  dt;  ds 


vdt, 


dx  = 


dy  = 


ds  = 


g- 


d^ 


-  tan  0  dO 
g 


g 


sec  0  dO 


(7) 

(8) 

(9) 

(10) 


Of  the  equations  (6),  (7),  (8),  (9),  (10)  there  is  but  one 
which  introduces  the  law  of  retardation  and  this  is  equation  (6), 
which,  as  it  connects  the  velocity  and  its  direction  is  evidently 
the  equation  of  the  hodograph  and  will  be  referred  to  as  such. 


The  integration  of  the  equation  of  the  hodograph,  permits 
the  immediate  reduction  to  quadratures  of  the  complete  solution 
of  the  principal  problem.  For,  if  v  may  be  expressed  as  a 
function  of  f>,  equations  (7),  (8),  (9)  and  (10)  will  on  the  substi- 
tution of  this  value  of  v,  express  dt,  dx,  dy  and  ds  in  terms  of 
^  and  d^*.     Equations  (7),  (8),  (9),   (10)  are  general,  and  the 


equation  of  the  hodograph  furnishes  the  data  for  solving  the 
direct  problem. 

The  Inverse  Problem. 
By  means  of  the  equation  of  the  hodograph  and  the  rela- 
tion that  it  affords  between  v  and  0,  it  is  clear  that  x,  y,  t  and 
s  may  each  be,  theoretically  at  least,  expressed  as  a  function  of 

0  alone;  and,  by  the  elimination  of  0  between  two  of  the  result- 
ing equations,  any  one  element  may  be  expressed  as  a  function 
of  any  other. 

If,  on  the  other  hand  experimental  firings  show  that  a 
relation  exists  between  any  two  of  these  elements,  this  relation, 
by  the  solution  of  an  inverse  problem,  makes  possible  the  de- 
tection of  the  law  of  retardation  and  the  other  elements. 
Suppose,  for  instance,  that/(x,y)  =  0  is  found  experimentally 
to  represent  the  equation  of  the  trajectory.     It  is  desired  to  find 

F  (v) 
\     and  the  other  elements. 

1  ta„.  =  g  =  y' 

„  ,    d  tan^  ^sec^^d<^^dV_    ,, 

^  dx  dx  dx'^     ^ 

3  f]    ~ '^      ^^^  ^^^'  ^  =  1+  y",  and,  hence,  from  (2) 

A  dx  ^  ,  dt'  y" 

^  Tr  =^  V  cos  0,  hence  ^~^  =  —  '^- 

dt  dx'  g 


Hence 


dt  =A       ^  dx 
\        or 


and  t—j\       —  dx 

^y  o     \  or 


Since  j.,  = ,  differentation  with  respect  to  t  gives 


df  g 

dx     d^x  ^  _  gy'"      dx 
dt   •  df^  y"^    •  dt 

dV 


2 

in  which  y'"  = 


dx^ 

Hence  ^  =  -  ^'''  =  -  ?L^^  cos  0  -  -     ^(^) 
"^^^""^    dt-'  2y"'  C     ''''^^^~       Csec^ 

and  from  this 

F  (v)  _        g       y"        /T+ y'^  y"'      ,       3,. 

"-^--2-y--\  ^      =-^^.v^cos^^ 


As  a  general  illustration,  suppose  that,  assuming  an  equation 

y  =  X  tan  c^  —  ^^^^ — 2    •  Q 
•^  '        2V  cosV 

to  represent  the  relation  between  y  and  x,  and  that  it  is  found 
by  experiment  that 

this  being  a  relation  obtained  from  the  plotting  of  the  observed 
values  of  V  and  X  f or  y  =  0,  (i.  e.  at  the  point  of  fall)  so  that 
V:sin2^^  ,,/       X\ 

gX            ^^      V"^'  C  )' 
Then,  representing  the  successive  derivatives  of  Q  with  respect 
to  X  by  Q',  Q",  Q'",  etc.,  we  have 
-  dy  _  ^„ g 


^'=lx-^-^-Fcosv(^Q  +  f-Q') 


y  = 


^(3Q'  +  3xQ"  +  |1q'" 


V^cos  . 

and  hence 

F  (v)  y'"      ,       3.       3  Q'  +  3  X  Q"  +  I  x'  Q'"        ,       ,  . 

-7s --  =  —  ^  .  V  cos  0  =  — ^ cr\T2^^ ^   .  V  cos  o 

C  2  g  2  V  cos  (p 

whence  the  law  of  retardation  is  known. 

Since  for  V==Oor-p=0,    Q   =   1,   (we    then    have    the 

conditions  in  vacuo,)  it  follows  that  'f'  (V, -tt)  must  vanish  with 


"Wc)  = 


V,  or  vanish  with  g  .     Hence  '/'  I  0,  ;;  )  =  0  and  'f'(Y,  0)  =-^  0. 

The  solution  of  the  inverse  problem  from  observations  in 
curved  and  high  angle  fire  is  particularly  useful  as  throwing 
light  on  the  laws  that  govern  the  flight  of  oblong  projectiles  at 
considerable  elevations,  which,  as  experimental  results  indicate, 
differ  considerably  from  the  laws  ordinarily  taken  as  a  basis; 
and  this  difference  seems  to  be  due  to  the  behaviour  of  the 
projectile  in  flight,  which,  while  systematic,  is  such  that  the 
resistance  is  not  a  function  of  the  velocity  alone. 

General  Character  of  the  Trajectory. 
Resuming  the  differential  equations  of  motion 
d  (vcos  d)  _  vF(v) 
"   d^  gC 


dt  =:  —  ^  .  sec  ^  d^ 
g 

dx  =  -  ^  d^^ 

dy  =  —  -^  tan  ^  d^^ 
g 

ds  = sec  ^  dO 

g 
1.     Let  it  be  required  to  trace  the  change  in  velocity  in  the 
trajectory  frorh  the  equation  of  the  hodograph,  the  latter  being 
placed  in  the  form 

dv  .  (F  (v)   ,     .    .\ 

a.  When  v  is  infinite,  ^tt.  is  infinite.     That  is,  the  rate  of 

d^' 

change  of  velocity  is  infinite  for  a  finite  change  in  6.  Hence  0 
has  a  limiting  value  whieh  designate  i.  This  is  in  the  ascend- 
ing branch. 

b.  Since,  so  long  as  0  is  positive,   Tn  is  positive,  it  follows 

that  V  diminishes  as  0  diminishes.  When  6  becomes  zero,  as  it 
does  at  the  summit,  we  have  at  that  point 

.  V  F(v) 


dy 

d^ 


o 

c.  Beyond  the  summit  0  is  negative  and  increases  in  abso- 
lute value,  and  -^  will  continue  positive  and  v  will  continue 

to  diminish  until  -rz  becomes   zero   and   changes   sign   which 
occurs  when 

and  at  this  point  a  minimum  velocity  occurs,  its  value  being 
such  that 

F(vJ  =-gCsin^.. 

d.  Beyond  this  point  -^  is  negative  and  v  increases  as  0 

diminishes,  until  0  approaches  — ^  in  value,  in  which  case  at 

the  limit 

F  (v)  ^  , 
gC        ^ 


or  F  (v)  =  g  C 

That  is  the  projectile  has  a  limiting  velocity  v'  such  that 
F  (vO  =  g  C 

This  is  interpreted  readily  as  indicating  that  the  velocity 
increases  until  the  acceleration  due  to  gravity  is  balanced  by 
the  retardation  due  to  air  resistance. 

2.  Let  it  be  required  to  trace  the  curvature  of  the  trajec- 
tory in  a  similar  manner. 

The  radius  of  curvature,  r,  of  the  trajectory  is  given  by 
the  equation 

V  d^^        v'^  ,, 

^=    y  =  -gCOsfl 

whence 

o 

r  = sec  0 

g 

a.  This  has  infinite  values  for  v  =  oo  and  for  sec  /?  =  oo; 

that  is,  for  inclinations  i  and ^-. 

Since  r  is  negative  the  curve  is  convex  upward.  Since  the 
radii  of  curvature  at  the  extremities,  (which  are  infinitely 
removed  from  the  summit)  are  infinite  it  follows  that  the  curve 
is  asymptotic  at  its  extremities. 

b.  For  a  minimum  radius  of  curvature 

dr  2  V  dv        .,        v'         ...       .,      ^ 

-.^  = -jn    sec  0 sec  0  tan  0=0 

dO  g  do  g 

whence 

dv  V   ,       ,, 

d^  =  --T^^^' 

but  5  =  vsec.(^^J    +sin.) 

Hence,  at  the  point  of  maximum  curvature, 

—  \  tan  0^.  =  V  sec  0^.  f  ^^  +  sin  0^.  ) 
2  V  gC  '/ 

or  -  sin  0^.  =  2  ^^  +  2  sin  ^, 

I.  '    n  2    F  (v) 

whence  sm  0^.  =  —  -^ ^'_,^- 

3     g  C 

Resume. 
The  point  of  minimum  velocity  occurs  when 

^  =  -  sin  .. 
gC 

and  that  of  maximum  curvature  where 

F  (v..)  2     .     . 


8 

Hence  they  are  both  in  the  descending  branch  of  the 
trajectory,  and  6  is  negative  in  both  cases. 

From  this  we  conclude  that  since  at  the  point  of  maximum 
curvature 

dv  V    ,  n 

^==— 2-tan^ 

and  since  tan  6  is  negative  and  hence  ^  positive,  the  velocity 

is  still  diminishing  and  hence  the  point  of  maximum  curvature 
lies  between  the  summit  and  point  of  minimum  velocity. 

Nomenclature  of  the  Trajectory. 

From  what  has  been  said,  it  is  clear  that  the  form  of  the 
trajectory  is  as  indicated  in  Figure  2. 

At  A,  V  =    00,  ^  =  i,   r  =   oo;  X  =  —  oo;    y  =  —  00. 

0;  y  =  0. 


At  B, 

V  = 

=  V, 

6  = 

^, 

r  = 

sec 

:^;x  =  0 

Ate, 

V   = 

=    Vo. 

6  = 

0, 

r  = 

g" 

;  X 

=  x„;  y  = 

At  D, 

V  ^ 

=     Vr, 

d  = 

0,, 

r 

= 

r„.  = 

— 

v' 

-^  sec  ^,. 

AtE, 

V   = 

=  v^, 

6      := 

0.r 

At  F, 

V   = 

=     Vo, 

,  ^  = 

-  (i)^ 

r 

Vo, 

'  seco^ 

g   -^ 

AtG, 

V   = 

=  v, 

6  = 

— 

2 

,r 

=  — 

oo: 

X;y  =  0 


That  portion  of  the  trajectory  above  or  near  the  level  BXof 
the  origin,  where  the  muzzle  of  the  gun  is  assumed  to  be,  is 
alone  considered  in  practice.     The  classification  of  trajectories  is 

based  on  the  values  of  — -^  and  ^  obtaining  in  this  portion  of 
the  trajectory. 

Thus  Direct  Fire  embraces  a  portion  of  the  total  trajectory 
which  portion  lies  in  the  vicinity  of  the  summit  so  that  both  ^ 
and  to  are  small  angles,  f  being  generally  not  above  15°. 

Curved  Fire,  which  includes  also  High  Angle  Fire,  em- 
braces all  fire  with  greater  elevations.  Curved  Fire  is  gener- 
ally limited  in  practice  to  velocities  lower  than  those  used  in 
Direct  Fire. 


As  limiting  cases  two  interesting  conditions  should  be  first 
considered.  These  are  rectilinear  motion  in  a  resisting  medium, 
taking  g  =  0;  and  motion  in  vacuo,  taking  C  =  oo. 


As  a  particular  application  the  case  in  which  F  (v)  =  A„V" 
is  considered,  and  is  of  great  value  in  connection  with  curved 
and  high  angle  fire. 


10 


In  Figure  3,  the  trajectory  as  used  in  practice  is  illustrated. 


Limiting  Cases. 

1.  g  =  0 

The  equation  of  the  hodograph  becomes 

vy 7)7  =  0  ;    or    0  =  <p  . 

d  (v  cos  0) 

Hence,  since  ds  =  vdt,  dx  =  cos  ^  ds;   dy  =  tan  <p  dx,  it 
follows  that 

X  =  s  cos  ^;  y  =  s  sin  c^ 
We  have  ds  =  vdt 

d^s_dv__F(v) 
^^^^  df     dt  C 

Hence 


t  -=  C 


./      V 


— ^v 

F  (vj 


—  vdv 

F  (vr 


C  (T  (v)  -  T  (V)) 


C  (S  (v)  —  S(V)). 


This  is  fire  at  an  angle  <p  with  the  horizontal.     In  the  case 
in  which  ^  =  0,  we  have  horizontal  fire  and  then 
t 
C 

X 

C 

y  =  0 


=  T  (v)  —  T  (V) 
=  S(v)  —  S  (V). 


2.  C  =  00. 

The  equation  of  the  hodograph  becomes 
d  (v  cos  d) 


& 


0 


or 


V  cos 


constant  =  V  cos  ip. 


11 

Hence  v  =  V  cos  f  sec  0 

Substituting  this  value  of  v  in  expressions  for  dx,  dy,  dt,  ds, 
we  obtain 

,,  V    cos    ^     ,     ,  n 

dt  =  —  ^  d  tan  0 

g 

,  V'COSV     ,    ,  n 

dx  = "  d  tan  0 

g 

dy  =  -^^d(tan^^), 

ds  =.  -  T^^l  sec^  0  dd 
g 

^_v:cosv^    ^ 
g 

in  which  {d),  or  simply  {d)  represents 

sec'  0  diO  =  \  (tan  d  sec  (^  +  log^  (tan  0  +  sec  0)) 

0 

Integrating  between  0  and  ^ 

V   cos   (p      ..  .  n. 

t  =  (tan  (p  ~  tan  0) 

g 

V"cos'>  .,  .       ,. 

x  = (tan  (f  —  tan  (i) 

g 

V'cosV  (,  ..  .  om 
y  =  -7. —  (tan"  <^  —  tan-  ry) 
•^  2g 

V'cosV  ..  ,        .... 

From  these  by  division 

t    _  sin  ^  /^  tan  0\ 

V  g     V        tan  <^/ 

X   _  sin2^/..  tan  6 \ 

Y'  ~    2g    I         tan^/ 

s   ^sin2^/   (<p) (d)    \ 

V'-'        2g    Vtan^       tan  o)' 

y   _  sin'V  (^  tan'  ^\ 

F  ~    2g    V^       tan>j 

Point  of  Fall. 
Substituting  y  =  0  in  this  last  equation  we  have  for  the 
point  of  fall 

tan'  d  =  tan'  ^ 

or  d  =  —  (^  =  —  CO 
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and  this  value  of  o  in  the  other  ratios,  gives 

T       2  sin  ^        ^      2  V  sin  ^ 

= ^,  or  T  =    ^ 

V  g    '  g 


X 

F  g 

^  ^  sin  2y^       {(f) 


sin  2  ^       V      V  sin  2  <^ 
^,  orX  = ^ 


g 


Also 


g: 


,  or  S  =  - 

tan  <p'  g 


=  1 


_  V^sin2^     JfJ_  _  ^ 


tan  ^ 


tan^ 


V 

w  ==  ^ 

Hence,  for  the  point  of  fall  the  following  ratios  are  unity 
in  vacuo 

gT 


From  these 


2  V  sin  ^ 

gX       ^ 
Fsin2^ 

Vo.     _ 

V  ~ 
tan  CO  _ 
tan^ 

X 


V  T  cos  (p'  ~ 

gV 
2Xtan^ 

Summit, 
On  substituting  ^  =  0,  in  the  expressions  fort,  x,  y,  s,  etc., 
it  is  seen  that 


t  = 


x„  = 


Yo  = 


s.  -- 


V  sin  <p  _    1 
g        ~    2 

_  V^  sin  2  cp  ^ 
2g 

V^  sin^  <P 


X 


2g 
V  cos'  (p 


~^gT  =^Xtan^ 


•  (^)  = 


g  '''2. 

v„  =  V  cos  cp 
The  following  forms  of  some  of  the  expressions  already 
given  are,  while  true  in  vacuo,  also  useful  in  practice  for  low 
velocities  and  heavy  projectiles: 
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Yo  =   g  g  T'  =  -g-  X  (tan  (p  +  tan  oj) , 

Vco  T  COS  w  .  V  T  cos  (p  =  X' 

x„  _  tan  OJ 

X       tan  (p  +  tan  (o 

y  =-    Y^  V^  —  x)  ^^^^  '^  +  ( '  1  )  ^^^  ^) 

Hence  for  approximate  determinations  of  co,  Vo, ,  y^,  from 
observed  values  of  V,  (p,  T  and  X,  we  may  proceed  to  find  oj 
from  the  relation 

s:  T' 

tan  OJ  =  ^ tan  ^ 

and  then  Va,  from 

_    X"     secjp  sec_a> 

X  tan  OJ 

X   = 

tan^  +  tan  oj 

Applications. 

1.      TO  FIRE  AT  A  POINT  NOT  ON  THE  SAME  LEVEL  AS  THE  ORIGIN. 

Given  V,  x,  y,  to  find  ^. 

V  ST  X  /  X  \ 

tan  -  =        =  tan  (p  —  ^  ^.T      .      =  tan  0    1 1  —  ^  )  • 
X  ^        2  V"  cos  (p  ^        X  ^ 

Let  (p^  be  the  quadrant  angle  of  departure  such  that  the 

trajectory  due  to  ^,  and  the  given  V  will  have  its  point  of  fall 

at  X,  0. 


Then 

^   _     V     ^m^Yx 

and 

X         sin  2  <p^ 
X         sin  2  ^ 

Hence 

tan.  =  tan  ^(1      '^  l^A 
^  V         sm  2  ^  / 

whence 

sin  2  c^  —  sin  2  c^,  =  2  cos'  cp  tan  - 

that  is 

sin  2  ^,  =  2  cos  p  ^.^           _^) 

COS  s 

_  2  cos  ^  sin  {(p  —  0 

cos  c 
Now,  from  trigonometry 

2  cos  c^  sin  (cp  —  s)  =  sin  (2  ^  —  .)  —  sin 
Hence 

sin  (2(p  —  0  -=  cos  .  ( tan  ^  +  sin  2  ^ J . 
From  which  2  <^  —  -~  and  <p  may  be  found. 
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2.      THE  PRINCIPLE  OF  THE  RIGIDITY  OF  THE  TRAJECTORY. 

This    * 'principle"   assumes   that  (p  —  e  =  ^, .     But  it  has 
been  shown  that 

.    ^  2  cos  ^  sin  {f  —  -) 

cos  ^ 
Distinguishing  ^  —  ^  as  ^'  we  have,  since  ^  =  ^'  +  ^, 

2  cos  (eg'  +  0  sin  y'  _  2  sin  ^'  cos  y'  cos  £ —2  sin'  ^'  sin  ^ 

cos  £ 


sin  2  ^^ 


cos  - 

=  sin  2  9?'  (1  —  tan  ^p'  tan  0 . 
The  assumption  that  ^'  -=  f^  —  Ms  productive  of  no  material 
error  so  long  as  tan  f  tan  -^  is  insignificant  in  comparison  with 
unity.  This  is  generally  true  in  direct  fire,  but  in  curved  or 
high  angle  fire  the  method  of  Application  1,  should  be  used  as 
far  as  practicable. 

3.      MEAN  VALUES. 

(a).     The  Mean  Ordinate,  h,  of  the  Trajectory. 


I    yds         (    y  sec 
P  ds  f  sed^  0  dO 


but 

y 

Hence 
Therefore 


V^'sinV 


F(^ 


tan"  ^\        T  V  sin"  ^ 

-^^j;andy„=— 2^- 


tan  (f 


Yo 


=  1 


tan'^i 
tan'  <p 


^       /^  (1  ~  ^an'^  0  cot^  <p)  see'  d  dO 


y.. 


P 


1- 


sec^  0  dO 

1  /  sec'  (p 
4  H^)tan^ 


cot'  (p     j 


Taking 


yo 


X  tan  (p,  it  follows  that 


h 
X 


If.  1    sec'  ^  ,    1       ,        \ 


Taking 


h 

T 


{9) 
1    .  1    /sec'  ^  .     \ 

y^  =  — ^ — ",  it  follows  that 
2  g    ' 

sin  2  ^  /,  1   sec' ^ 

^  \^tan  (p  —  ^ 


4g 


{<P) 


f  cot^). 
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The  mean  height  is  useful  in  correcting  for  altitude  and  is 
here  given  in  terms  of  X  and  <f,  and  of  V  and  c^.  The  latter 
gives  too  large  a  value  for  h  but  does  very  well  (when  X  is  not 
known)  as  a  first  approximation. 

(b).     Mean  values  of  cos  0  and  of  sec  (f. 

1  '  cos  d  sec'  0  &)     _     tan  c^ 

1  (cos  ^V),„  =  •' °         (^^  J^ 

I  '  sec  0  sec''  0  &)  _  3  tan  <p  +  tan''  ^ 

2  (sec  ^y),„  =  •' °       ^^^-  ~  3  (^)    -- 

3  («-e«-i).'  =  sec^ 0',„  =   '    \ ,,  =  3tan<.  +  tanV  ^ ^  + .  t^„. 

(cos  0),,  cos^,„  stance 

or  l  +  tan"' 6'',,,  =  l  +  itan"'c^ 

Hence  tan  /^',„  =  ±         ' 

I    o 

(^'n,  is  thus  an  angle  for  which  (cos  /^),„  and  (sec  ^^),„  will  consist. 

It  is  the  best  determination  for  general  use.     The  double  sign 

does  not  signify  as  ^',„  is  used  to  obtain  only  cos  ^^',„  and  sec  6/',,, 

which  are  even  functions  of  tan  o. 

4.      DROP  DUE  TO   GRAVITY. 

(a).    Jump. 

A  projectile  emerging  from  the  muzzle  of  a  piece  drops, 
due  to  the  action  of  gravity,  from  the  straight  line  y  =  x  tan  <^, 
on  which  it  departed,  a  vertical  distance  given  by  the  equation 
of  the  trajectory  in  vacuo;  thus: 

2f  X^ 

y  =  X  tan  ^  —  ^'  ^r->   — •' 
V"  COS'  <p 

2f  X' 

Hence  the  drop  is  v  ^y.,^ — .  -  feet. 

V'  COS"  ip 

The  gun  being  laid  by  bore-sights  on  a  mark,  the  height  at 
which  the  projectile  actually  strikes  must  be  increased  by  the 
drop  in  order  to  give  a  point  in  prolongation  of  the  line  of 
departure. 

This  gives  the  means  of  determining  the  amount  and  direc- 
tion of  the  jump. 

The  drop  in  inches  is  ^^fr^^> — 
V  COS"  f 

in  which  as  before  x  is  the  distance  in  feet  to  the  screen  on 

which  the  jump  is  measured. 

(b)     Danger  Space. 

Using  the  expression 

y  -=  ^  (X  —  x)  (x  tan  oj-\-  (X — x)  tan  if) 
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and  representing  X  —  x,  the  danger  space,  by  D,  there  results 

y  =  '^^.  D.  ((X  — D)tanoj  +  Dtan^) 

or  y  X'  =  D  (X  —  D)'  tan  co  +  W  (X  —  D)  tan  <p 

whence  dividing  by  X'^  tan  co 


X  X        ^^  '  ^^  V  X  /     '   ^  V  X 

tanj^ 
tan  oj 


in  which  P  =  1  — 


A  table  of  values  of  ^-^ —  may  be  computed  with  P  and^ 
as  arguments,  and  ^  taken  out  by  interpolation  using  ^ — 
and  P  as  arguments.     From  this  table  may  be  computed 


ycotw 


-p  •        -o      J  y  cot  w 

for  a  given  P  and  ^ — ^ — . 

Ordinarily,  however,  the  last  term  is  insignificant  and  the 
second  comparatively  so,  and  the  following  method  may  be  used: 
Take  D  =  y  cot  ^  as  a  first  approximation  and  substitute 
this  value  in  the  last  two  terms,  thus 

ycoto^  ___D  _  .p  , -,>,  y'  cot" OJ      p  y' cof  qj 

X      "X      ^^^^^     r     ^^     r 
/p+i  \  p 

D  =  y  cot  w  +  \x)^'  ^^^'  ^  —  X^  ^'  ^^^'  ^  =  X  —  X 

in  which  P  =  1  -^^ 

tan^ 

This  formala  may  be  used  to  find  x,  when  y  is  given,  on  a 
given  trajectory  especially  in  curved  and  high  angle  fire,  and 
is  useful  for  that  purpose.  Ordinarily,  for  danger-space  prob- 
lems P  may  be  replaced  by  0  giving 

^  ,      ,  y'  cot'  OJ  /^    ,    y  cot  ^  \ 

D  =  y  cot  "^-^      Y —  ^y  ^otoj[l  +  ^  ^      j. 

5.      MAXIMUM  AND  MINIMUM  VALUES. 

(a).     Range,  as  a  function  of  Elevation. 
X  =  V  sin  2  (p 

T~-  =-   _ .  2  COS  2  c^  =  0  for  a  max. 


Hence  cos  2  ^  =  0;  2  ^  =   ^   ;  <P  =  ~~r  =  ^5 
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Hence,  in  vacuo,  45°  is  the  angle  of  greatest  range.  This 
is  also  nearly  true  in  curved  fire,  the  value  of  <p  being  some- 
what under  45°. 

(b).     Time  of  Flight,  as  a  function  of  Elevation. 
^  _    2  V  sin  ^ 

g: 

dT       2V 

T     =  —  COS  (f  ^  y) 
d^        g         "^ 

Hence  cos  ^  =  0  and  ^  =  -^  =  90°. 

(c).     Minimum  velocity  in  the  trajectory  due  to  V,  <p. 
V  =  V  cos  ^  sec  0 

-To   =  V  COS  ^  sec  0  tan  0  =  {)  for  min. 

Hence  sec  6  tan  ^^  =  0  or  ^  =  0. 

Hence  minimum  velocity  is  at  the  summit.  This  is  nearly 
true  in  curved  and  high  angle  fire  with  low  velocities,  but  is 
not  even  approximately  true  in  direct  fire. 

(d).     Maximum  Curvature,  as  a  function  of  0. 

ATcr   X.         v'  ^  V"'  V  cos"  f  sec'  0 

We  have  —  =  —  g  cos  0  or  y  =  — ^  =  — ^h 

r  g  cos  0  g  cos  0 

TT  V'  cos'  <p  3  a 

Hence  r  =  — ^.  sec  o 

§: 
and 

d  r  V  cos'  (p     ^        ^  a  .       n       n  ^ 

j-^  = —  .  3  sec  o  tan  0  =  0  lor  mm.  r  or  max.  curva- 

d^  g 

ture,  which  gives  ^  =  0. 

or  the  maximum  curvature  occurs  at  the  summit. 

This  is  not  absolutely  true  in  air,  but  is  nearly  so  in  curved 
and  high  angle  fire. 

These  maxima  and  minima  serve  as  points  of  departure  in 
seeking  those  in  air,  rendering  fewer  approximations  necessary. 

DIRECT    FIRE. 

This  class  of  fire  pertains  to  small  angles  of  departure, 
generally  not  above  15°,  though  it  may  be  used  for  greater 
elevations  with  the  introduction  of  a  suitable  factor  in  the 
ballistic  coefficient. 

The  method  here  used  is  that  of  Siacci,  in  which  is  used 
the  preudo-velocity  which  is  designated  u,  and  which  fulfils 
the  condition 

u  cos  <p  ^  V  cos  0 
or  u  =  V  cos  d  sec  ^ 
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The  integration  of  the  hodograph  is  accomplished  by  the 
use  of  an  integrating  factor,  ^5,  such  that 

F  (v)  =  /?.  F  (u)  cos'  <p  sec  6, 

with  the  practical  certainty  that,  for  direct  fire,  the  value  of  /? 
will  never  differ  greatly  from  unity.  In  any  case,  however,  we 
may  write  the  equation  just  assumed. 


C 

in  which  C, 


C 
C 


F  (u)  cos'  (f  sec  d 


FJu)^ 


cos^  (f  sec 


We  will,  however,  use  C  instead  of  C,  with  the  understand- 
ing that  where  Siacci's  method  is  used,  C  is  to  be  divided  by  /?. 


IS  a 


Figure  3  shows  the  character  of  u.  It  is  thus  seen  that  u 
IS  a  velocity  parallel  to  the  muzzle  velocity  and  having  the  same 
horizontal  component  as  the  actual  velocity  v. 

The  equation  of  the  hodograph  becomes,  under  these 
conditions 

^  (u  cos  ^)    ^  u  COS  f  sec  6.  F  (u) .  cos'  (f  sec  0 


or 


d^ 

du       _  cos'  f  sec'  6  dd        cos'  ^      ,  ^      n 
u  F  (u)  ^C  ~   TC~  •  ^  t^^  ^ 


or,  integrating  between  6  and  f 

tan^— tan^  =  ^^-^-  /     ^ 
2  cos  <p  i  u 


—  2gdu 
F  (u) 


£v(^^^)-^(^)) 


2  cos'  iP 
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The  expressions  for  dt,  dx,  dy,  may  be  written 

,,  u  cos  ^     -,  ,       ^  ^  du 

dt  -= —.  d  tan  0  =  —  C  sec  ^  -^-/^ 

g  "^  F  (u) 

,  -,,  Cu  du 

dx  =  u  cos^  dt  =  —  p  (^)^ 

dy  =  ucos^tan^dt  =  — -^^~^(tan^  — 2^^  (I(u)— I(V))j 
and  from  these 


t  =-Csec^/  — -^~V    =Csec^     T  (u)  -  T  (V) 


I  (V))dS  (u). 

2^  (a(u)-A(V)-I(V)(S  (u)-S  (V))) 
in  which 

A(u)  =/l(u)  d  S  (u) 

X 

Noting  that  S  (u)  —  S  (V)  =   -p-  and  dividing  both  members 

of  the  equation  for  y  by  x, 

y       ^  C       /A(u)-A(V)       J  ,^,A 

tan  ^  =   "^  -  =  tan  f  —x ^  I  ^7^^ — ^.^~  —  I  (V)  ) 

X  "^       2cosV  ^S(u)  — S(V)  / 

The  angle  -  as  thus  defined  is  the  position  angle.  In  Part 
I,  its  absolute  value  alone  was  considered,  but  here  it  is  of  the 
same  sign  as  y. 

The  S,  A,  I  and  T  functions  of  u  constitute  the  principal 
ballistic  functions  and  are  in  general  use  in  all  the  principal 
countries  in  the  world.  From  a  table  of  these  functions  certain 
combinations  were  derived  by  Major  Braccialini  and  tabulated. 
These  functions  materially  abridge  the  labor  of  computation, 

and  are  tabulated  with  V  and  -p-  as  arguments.     Of  course  the 

values  of  all  these  functions  will  depend  on  the  law  of  resis- 
tance adopted.  That  in  use  in  this  country  and  forming  the 
basis  of  Ingalls'  Tables  I  and  II  is  due  to  General  Mayevski,  is 
based  on  the  Krupp  firings  at  Meppen  and  was  extended  by 
Colonel  Zabudski  to  include  velocities  up  to  about  3600  f-s.* 

*  See  Appendix  III. 
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Braccialini's  secondary  functions  are  as  follows: 
t'  =-^cos^  =  T  (u)  -T  (v). 
,       2  cos'  f  (tan  ^  —  tan  0) 


^  =I(u)— I(V) 

_  2  cos'  f  (tan  ^  —  tan  ^)    _  A(u)  —  A(V)  _  ^  .^j. 
^  C  S(u)-S(V)       ^^^ 

,    _  2  cos  <p  (tan  e  —  tan  6)  _    ,  _ 

D   —  y-^  —   3,  Si 

z  =  ^  =  S  (u)  -  S  (V) 

,  _  2  cos'  (f  (tan  (f  —  tan  e)   ^  _b^ 
X  z 

.,  _  tan  ^  —  tan  e   _  Jc^ 
tan  e  —  tan  0  a 

The  type-symbols  of  these  functions  are  not  in  all  cases  the 
same  as  those  used  by  Braccialini,  but  are  those  used  by  Colonel 
Ingalls  and  are  here  given  so  that  there  will  be  no  confusion 
in  using  Tables  I  and  II,  Artillery  Circular  M. 

Using  the  secondary  functions,  working  equations  take  the 
form 

a'C  ^  /,  a/C    \ 

tan  (p  —  ^ 2 —  ^  tan  ^  ( 1 -. — ^ —  J 

2  cos  <p  ^         sm  2  ^  / 

t  =  Ct'  sec  <p 

y        .  a  C  ,  /-,  a  C     \ 

—  =  tan  (p  —  ^ Y  —  =  tan  ^  ( 1  —     .    ^ —  ) 

X  2  cos  (f  \  sm  2  v'  ^ 


tan 


tan  e 


tan  <p  —  tan  £    _  ,  , 


tan  £  —  tan  0 

sin  2  (ff    tan  £  \   _    , 

X      ^        tan  cp) 

X  =  Cz 

The  function  a'  has  another  use,  shown  by  Colonel  Ingalls. 

It  is  as  follows: 

For  the  summit,  6^  =  o;  hence,  at  the  summit, 

,  _  sin  2  (p 
a   -       ^ 

It  is  also  true  that,  at  the  point  of  fall,  where  y  =  o, 
_  sin  2  cp 

^  ~  ~~^ 

Affecting  the  letters  referring  to  the  summit  with  a  zero 

subscript,  and  using  the  corresponding  capital  letters  for  the 
point  of  fall,  these  two  conditions  are  expressed  as  follows: 
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The  condition  a,/  =  A  at  once  determines  Sid  when  A  is 
known,  and  Col.  Ingalls  used  this  relation  in  determining  a 
function  ao"  for  the  determination  of  the  maximum  ordinate. 
Thus 

^'-  =  tan  <p  (l  -^"  )  =  tan  <p  ^^^^  =  ^'^  tan  <p 
Xq  ^         a  ^  ao  ao 

whence 

yp      ^  ^0     jOq^  ^  Op  Zq  =  o  " 

Ctans^        C    •    ao'  a^' 

So  that    yo  =  ao"  C  tan  cp. 

Thus,  in  using  the  tables,  ao'  is  at  once  known  when  A  is 
known  and  all  formulas  relating  to  the  summit  may  be  used  on 
taking  out  the  necessary  quantities  by  interpolation  using  ao'  as 
argument. 

The  calculation  of  the  functions  a,  a',  a",  b,  b',  c',  t'  is 
undertaken  with  assumed  values  of  z  at  regular  intervals  and 
with  assumed  values  of  V  at  intervals  allowing  satisfactory 
interpolation. 

Thus  with  assumed  values  of  V  and  z,  u  is  at  once  found 
by  means  of  Table  I  from  the  relation 

S(u)  =  z  +  S  (V). 

The  S,  A,  I  and  T  functions  of  u  and  V  are  to  be  found  in 
Table  I  and  the  excess  of  the  functions  of  u  over  the  corres- 
ponding functions  of  V  are  found.     Representing  these  excesses 


by  As,  Aa.  Aj. 

At,  so  that 

As    =  S  (u)  —  S  (V) 
Aa    =   a  (u)  -  a  (V) 
Ai     =   I  (u)  —  I  (V) 
At    =   T  (u)  —  T  (V) 

it  follows  that, 

a  =  —  Aa   —I  (V) 

Zi 

a'  =  A, 

t'  =  At 

The  remaining  functions  are  readily  found  when  a', 

a  and 

z  are  known,  thus 

b  =  a'  —  a 

b'^    b 
a 

z 

^    ~   a' 
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Point  of  Fall. 
The   equations  for  the  point  of  fall,  where  y  =  o,  become 

X  =  cz 

T  =  C  T'  sec  ^ 

sin  2  ^  =  A  C 

tan  (jf)      -r>r 
=  ^ 

tan  (p 
Yu>    =  u<o    cos  (p  sec  CO 

Su7nmit. 
Using  ao'  =  A,  take  out  Zo,  slq"  and  any  other  desired  func- 
tions for  the  summit,  then 

Xq  =  Czq 

yo  =  ao"  C  tan  <p 

to  =  Cto'  sec  (p 

etc. 

Any  Point  on  the  Trajectory, 
tan  ^  =  ^^^^  (a- a') 


A 

tan  £  =  -^ 


y  =  toM^  A 

X  A      V 


The  other  elements  may  be  found  in  the  same  manner  as 
for  the  point  of  fall  using  the  proper  value  of  z  as  argument. 
It  is  a  property  of  pseudo-velocity  that  for  a  given  V  and  Z  it 
is  independent  of  the  angle  of  departure  (p  or  of  the  height  y. 
This  is  true  because 

S  (u)  =  z  +  S  (V) 
showing  that  u  is  dependent  on  z  and  V  alone,  and  since  the 
secondary  functions  are  dependent  only  on  Z  and  V,  they  are 
all  independent  of  ^  or  y. 

Represent  by  ^^  the  value  of  (p  which  would  for  a  given 
V  and  z  cause  the  point  of  fall  to  occur  at  x,  o,  the  foot  of  the 
ordinate  y ,  whose  abscissa  is  x.  In  that  case  e  becomes  zero  and 
hence  sin  2  ^^  =  a  C 

It  is  also  true  that 

a  C  =  2  cos'  (p  (tan  (p  —  tan  £)  =  sin  2  ^  ( 1 , j 

Hence 

sin  2  (p^  =  sin  2  ^  —  2  cos'^  tan  £ 
This  is  exactly  the  equation  already  found  for  conditions 
in  vacuo,  and  hence 

2  cos  (p  sin  (cp  —  e) 


sin  2  (p^ 

cos  s 


and 
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sin  {2f  —  e)  =  cos  e  (tan  £  +  sin  2  ^J  =  cos  £  (a  C  +  tan  e) 
which  may  be  solved  for  f,  when  V,  x,  y  and  C  are  given. 

Since  the  equations  are  identical  in  form  with  those  for 
conditions  in  vacuo,  the  error  in  elevation  due  to  the  assump- 
tion of  the  principle  of  rigidity  of  the  trajectory  may  be  found 
from  the  relation 

sin  2  ^^  =  sin  2  ^'  (1  —  tan  f  tan  e) . 

Forms  Analogous  to  those  in  Vacuo. 
The   equations    already   found   may  be  placed  in   forms 
analogous  to  those  in  vacuo;  thus 


cr  x^  /  ST  X  Q         \ 

y  =  X  tan  <p  —  ,7^^ ~  %  =  ^  tan  ^  1 1  — ^yr  -  ^d—  ) 

2V  cos^   •'  V        V  sm  2  ^  / 

md  =  tan  ^  —  -;ry;,~ — ^ —  .  q^  =  tan  c^    1  —  f^f-^-^  ) 
^       V  cosV  '    V        V'sm2^/ 


t=    ^ 


Vcos^* 

V  COS  ^  =  V  cos  (f.  q,,. 
in  which 

V^  sin  2  ^         a  C  a  V^        F  c 

^  g  X  sm  2  ^  g  z  g 

_  V'  sin  2  ^       a  C     _  a'JVi 
^'  2gx      *  sin2>       2gZ 

V  cos  (P        n  Vt' 
q,  = ^^  .  ct  sec  f  =  — 2" 

V  cos  d  sec  <p  u 

q.  = V —  =  V 

The  form  q,  is  exceedingly  useful  in  connection  with  experi- 
mental firings  to  determine  c.  It  varies  less  rapidly  with  V 
than  does  c'  and  interpolation  is  facilitated  by  its  use. 

All  these  functions  are  especially  useful  in  the  case  in 
which  F  (v)  =  A„  v",  as  will  now  be  shown. 

CASE   OF   F   (v)   =  A„  v°. 
Representing  by  m  the  product  A„  V"  ~  ^  z  and  noting  that 

S  (u)  -  S  (V)  =  ^^^^^li^-ll   )  =  z 
T  (u)  -  T  (V)  =  („„i}a,.V--^((¥)  "  -  0 
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it  follows  that 

q»=  — .    ((1+  (n  — 2)  m)^'^2--l') 
n  m      V  / 

-  1 

q,  =  (1+  (n  —  2)  m)n-2 

The  principal  entire  values  of  n  are  2,  3,  4;  and  for  these 
the  following  forms  are  taken  by  the  functions. 


n 

2 

3 

4 

m 

A,  z 

A,  Vz 

A,V^z 

Qy 

e^"  —  1  -  2m 

l  +  f  m+^  m' 

1  +  1  m 

q^ 

e'^-^-l 
2m 

1  +  m  + 1  m' 

1  +  m 

Qt 

e"^  —  1 
m 

1  +  1  m 

(l  +  2m)^^-l 
3m 

Q—  m 

1 

1 

Qv 

1  +  m 

(l+2m)^ 

From  the  expression  for  q^  it  is  seen  that  the  trajectory  is 
a  curve  whose  degree  is  — — ^,  and  hence  the  degree  of  the 
trajectory  for  the  following  values  of  n  is  as  indicated. 


n 

00 

0 
1 
2 
3 

4 

8 
3 


degree 

2 

* 

00 

4 
3 
5 
6 

7 
8 


F(v) 
0 

A., 
A,  V 
A.  v'^ 
A^v^ 
A,v* 


Condition 
Vacuum 

Const,  resistance 
Linear 
Quadratic 
Cubic 
Biquadratic 


2m 


m 


m 


A.v" 


It  will  be  noted  that,  for  n  =  2  and  below,  the  trajectory  is 
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transcendental.     For  n  =  2   the   familiar  limit    (l  +  p)^p  on  p 
approaching  zero,  occurs,  and  is  in  the  form 

(1+  (n  — 2)m)n-^  I      .  _  =  the  Napierian  base  e. 

In  the  case  n  =  1  and  n  =  0  illusory  forms  occur;    but,  on 
evaluation,  show  logarithmic  forms.     Thus: 

n  0  1 

m  '  Ao  A,  V 

q,      -  4((l-2m)log.(l-2m)-2m)       -  ^(l+  lo^^-^l-"^)\ 


q^      -  i  log,  (1  —  2  m) 
1 


m 
1— m 


q,      -^((1  — 2m)i  — 1  )  —  mlog,  (1  — m) 

q^  (1  —  2m)^  1  —  m 

Expansion  in  series,  for  the  general  exponent  n,  gives  the 
following: 

4  4  o       4 

q,  =  l+  r^   m  +  r  .-  (4  —  n)  m"  +  ■  _    (4  —  n)  (6  —  2n)  m^  +  etc. 

II.        I±  lA 

2  2  2 

qe  =  1  +  r^  m  +  j       (4  —  n)  m'  +  jx  (4  —  n)  (6  —  2n)  m'  +  etc. 

q,  =  l+  .      m  +  .-g-  (3  — n)  m'  4-  .  ,    (3  —  n)  (5  —  2n)  m'  +  etc. 

q,  =  l— pp  m  —  r^  (1 — n)  m"  —  r^  (1  —  n)  (3  —  2n)  m^  +  etc. 

These  series  will  have  a  finite  number  of  terms  for  values 
of  n  as  follows,  representing  the  number  of  terms  in  the  series 
by  r: 

^  ,  2  r      .  2r— 1    .  2  r  —  3 

Forq.andqe  ,n  =  ^_^;  f or  q„  n  =  — :^^;  forq,,  n  =    ;^^^^  • 

THE  CONSTRUCTION  OF  TABLES  FOR  THE  CASE  F  (v)  =  A„  V". 

With  m  and  n  as  arguments  the  values  of  the  q-factors 
or  their  logarithms  may  be  readily  calculated,  and  tabulated, 
together  with  such  additional  factors  as  may  be  derivable  from 
them,  using  the  formulas  already  deduced. 

There  is  one  case  remarkable  on  account  of  the  fact  that  V 
is  not  involved  in  m.  This  is  the  case  in  which  n  =  2,  and, 
as  a  consequence, 

m  =  A.,  z 

This  is  useful  in  curved  fire  as  it  is  nearly  true,  for  the  low 
velocities  used,  that  the  quadratic  law  obtains.     Knowing  A.  we 

may  use  -r-  or  z  as  the  sole  argument  of  the  table. 

A2 
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Quadratic  Law  of  Resistance. 


Placing 


Qi  =  g  q^;  q2  = 1;  qs  =  k^;  q.  ==  ". 

qa  q2  1      1 

qs  =  z  q^;  q^  =  ~^'  ;  q,  =    »   .      ^^g,  q^ 
q*  -^  -^2 


We  then  have 

y  =  X  tan  ^ 


X"  q, 


2V^cosV 

X  qi  q^+l 


tan  d  =  tan  ^  —  ,,.,       ^      .      « 
^       V"  cos  (p         2 


t 

V 

2F 

X 

qi 

V  cos  ^ 

V  cos  (p  sec 

2q3 

q4 

cos  (p  sin  (^ 

-0 

2  F 

C  cos  £ 

COS  /^  sin  (s  - 

-^) 

q5 

C  cos  £  ~   ^   "^^ 


Zo  =  Q 


Yo      ^   Zq/-,      Zq     qi  (Zq)  a    p. 

X  tan  ^        Z  V^       Z  •       Q,     /     ^' 
in  which  q^  (zj  is  the  value  of  q^  for  z  =  z„  and  Qi  refers  to  the 
point  of  fall, 

POINT  OF  FALL. 

At  the  point  of  fall  £  becomes  zero  and 

V  sin  2  ^ 


X 

tan  00 

tan  ^ 

V  sin  ip 

T 

V 


Qi 
Q2 
Q3 


V^  sin  2  ^  _  ^ 

C  "^^ 

v^to  sin  2  <^  _  p^ 

C  ~  ^' 

Zq    ^    V3i7 

X  tan  9      ^' 
These  together  with  the  equations  for  the  general  point,  x, 
y,  of  the  trajectory  will  give  a  complete  solution  for  this  law. 
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ANY  POINT  X,  y,  OF  THE  TRAJECTORY. 

Since,  as  already   shown  for  Siacci's  method,  the   factors 

depend  on-^  and  V  alone,  it  follows  that  the  solution  for  any 

point  may  be  readily  accomplished  by  the  use  of  the  equations 
already  given  and  by  the  equations 

sin  {2<p  —  s)=  cos  e  (tan  ^  +  sin  2  ^ J 

.    o  2  cos  (f  sin  {(f  —  e)         X 

sm  2  ^x  = =  -77-2 .  qi 

cos  s  V^ 

From  these  and  the  preceding  are  found 


y        ,  /^       sm  2  ^^  \ 

-^  =  tan  ^  ( 1 ^~K 

X  \  sm  2^  / 

tan  ^  =  tan  ^  fl  -  (q.  +  l)  ^I?A^)  . 
V  sm  2  ^  / 


V      sin  2  (f, 
t  cos  <p  = 


2   •         qa 

V  cos  if  sec  t 


1  0 


V  = 

q4 

The  Integrating  Factor. 
It  must  be  remembered  that  the  relation 

F  (v)  =  /?  F  (u)  cos-  f  sec  0 
becomes  for         F  (v)  =  A„  v", 

A„  v"  =  /?  A„  u"  cos-  f  sec  0 
so  that 

/?  =  ( — j  cos  0  sec^  f  =  cos"  "^  ^  sec" 

For  curved  fire  it  is  necessary  to  ascertain  and  use  the 
value  of  /?.     Using  n  =-  2  we  find 

/?  =  sec  6^ 

Its  mean  value  is  readily  found  over  the  range  {0),     Thus 

V 

cos  0  sec**^  0  dd 

oj tan  <p  +  tan  ^ 

Finding  also  the  mean  value  of  P  in  order  to  determine  a 
mean  mean  value  for  /?,  thus 

j       sec'  d  d^ 
J  —CO 3  (tan  <p  +  tan  co)  +  (tan^  ^  +  tan^  co) 

{<p)  +  {co)        ~  3((^)  +  (o;)) 

we  have 


.  / 


Hence 
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3  (tan  <p  +  tan  <^)  +  tan^  ^  +  tan^  ^ 

3  (tan  f  +  tan  oS) 
1  +  .T  (tan^  ^  —  tan  ^  tan  u)  +  tan^  oS) 

tan  <^  ,  tan^  o)  \ 


/5  =  Vl-^-tan'-(l-i^> 


tan  ^  '  tan^  ^  / 
A  fair  mean  value  of  /?  may  be  gotten  by   assuming  for 

curved  fire  the  mean  ratio  , =  '6/5,  which  gives  practically 

?  =  v/X^^ftan'V 
This  applies  with  great  accuracy  to  the  tables  of  Curved 
Fire,  giving  results  almost  identical  with  those  obtained  by 
quadratures.  It  applies  to  all  elevations  up  to  about  60°  and 
performs  its  function  of  producing  the  same  range  with  the 
same  V,  ^,  and  projectile  when  computation  is  made  by  the  two 
methods.     It  has  no  other  function.  «. 

Differential  Formulas:    General. 
It  is  often  desirable  to  find  what  effect  a  change  in  V,  (f,  or 
C  will  produce  on  the  range  to  the  point  of  fall. 
The  relation 

^^Mu)^AiV)__i(V) 

or  (A  +  I  (V))  Z  =  A  (u)  -  A  (V) 

gives  on  differentiating  with  respect  to  A,  Z,  and  S(V), 

I  (u)  d  Z  +  A'  d  S  (V)  =  Z  d  A  +  A  d  Z  +  I  (V)  d  Z  +  Z  d  I  (V) 

or  (A'— A)  d  Z  +  A'  d  S  (V)  =  Z  (d  I  (V)  +  d  A) 

But  dI(V)  =  ~f-dS(V)  =-^|^^^ 


/I  Q 

r\T\ 

VdV 

CI  o  V  V  ;    — 

F(V) 

A'- 

-A  =  B 

ence 

BdZ 

A'  V^ 
F(V) 

d  V 
•   V 

nence 
dZ 

A'  V^ 

dV 

Z 

ZBF(V) 

•  V 

A'V^  — 2 

:gz 

Z  B  F  (V)     • 
Now  Z  =  -p  and  hence      ^  ^ 

also 


z(d  A 

2g 

F(V)  • 

dV 
V 

dA 

2g 

^      B 

BF(V)    • 

dV 

dA 

V      ' 

B 

dZ 

dX 

dC 

)■ 

_dVl 
V 


d  A          A   /d  sin  2  <^ 
A"         B  ^"sin2^ 

dC 
C 

sin  2  <p 
C 
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d  A    _   A 
B  B 

since  A  == 

Hence 

_dX   _dC   _  A'V^  —  2gZ       d^V^  A_yd  sin  2  ^  _  d  C\ 

X  C  Z  B  F  (V)      •      V  B  V  sin  2  c.  C    / 

or 

dX^  _   (    _P^\dC        A      dsin  2w  A'  Y'  —  2gZ       dV 

X  V^       B/C  B    •    sin27  Z  B  F  (V)       *     V 

1    .         A       tan  (f         J  A'      ^    ,  tan  o       , 

and  since   -r»  =  j. ,  and  ^h"  =  1  +  ^ —    and 

B       tana>  B  tan^ 

-r»        A  T^/       sin  2  c.     tan  oj         .      ^,        tan  a; 

B  =  A  B    =  —y,—^  .  ^ ;   also  B    = 

C  tan  (f  tan  (p 

the  equation  becomes 

d  X  _  /.  _taiLf \  dC 
X        V        tan^o/"C' 
tan^^     d  sin  2^ 
tan  CO  '    sin  2  ^ 

C      V^     (.    ,  tana._     2gX  Uan  s^     d^ 
X  F  (V)   V         tan  <p       V^'sin2^/tan  co  '    V 
which  may  be  used  directly  with  a  range-table  if  tabulated 

values  of  t^  /ttn     are  at  hand,  or  if  the  value  has  been  found 

for  the  particular  range-table. 

The   method  of  finding  range  changes   by  increments   is 
given  in  Part  I. 

Special  Differential  Formulas  for  Range  Changes. 
Taking  the  development 

^         V  sin  2  ^       ^    ,    o        ,     , 
Q,  = gir^    =  1  +  t  m  +  etc 

X 

restricting  it  to  two  terms  and   noting  that  m  =  A„  V"  ~  ^ 

it  becomes 

gX         1  +  ^^"^         c 
Differentiating,  it  is  found  that 

V-sin2^/2dV      dsin 2^     d  Xx   _  _X  /  .dV     dX 

gX      V  V    +    sin2^        X  J   "^^'^^       C    V^""    ^^  V  "^  X 
dC\      /T-sin2^      ,W,         ^.  dV^dX      dC\ 
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Representing  ^   ^  —  1  by  N,  and  collecting  terms 

(2+  (4  -  n)N)  ^  +  (N  +  1)  ^^'^  +N^  =  (2N  +  1)  ^ 
\  /    V  sm  2  ^  C  X 


or 


dX_  f^_   nN   \dV  ^    N  +  1    dsin2^  ^      N     dC 


/         nN 

V^     2N  +  1/ 


X  V"  2N  +  1/  V  2N  +  1*  sin2^  '  2N  +  1  C 
nt  is  approximately  2  for  long  ranges  with  muzzle  velocities 
above  1800.  For  field  guns  with  velocities  about  1700,  n  =  3^^. 
In  general,  for  a  range  table,  n  may  be  taken  as  k  + 1  R,  R  being 
the  range  in  yards;  and  k  and  1  constants  giving  proper  X 
changes  due  to  V  changes  at  long  and  short  ranges. 

In  these  formulas  it  is  necessary  to  know  only  V,  ^, 
and  X,  in  order  to  apply  them,  provided  n  is  also  known 
approximately. 

For  curved  fire  n  =  2  giving 
dX_2N  +  2     dV  ,    N  +  1      dsin2^_^      N        dC 


X       2N  +  1  •    V       2N  +  1  •    sin2^       2N  +  1  '     C 

For  field  Artillery,  V  between  1000  and  1700,  n  =  31,  hence 
dX  ^   N  +  4      dV       N  +  1      dsin2^  N         dC 

X       4N  +  2  •     V       2N  +  1  •    sm"2'^      2N  +  1-    C 

The  differential  formulas  apply  to  small  changes  in  range 
due  to  changes  in  the  different  elements.  They  are  useful  for 
quick  computations  using  only  a  range  table  and  a  table  of 
logarithms;  they  are  also  of  value  in  showing  what  will  result 
if  a  change  is  made  from  one  gun  to  another,  as,  for  instance, 
when  it  is  desired  to  find  what  change  is  to  be  made  in  muzzle 
velocity  when  the  gun  and  projectile  are  changed,  in  order  to 
acquire  the  same  range  with  the  same  angle  of  departure.  In 
that  case  (using  for  illustration  our  approximate  formulas  with 

properly  determined  n) ,  -^^-  and  —. — ^      become  zero  and  hence 

X  sm  2  ^ 

/2+(4  — n)  N  \dV  N      dC   _  . 

V        2N  +  1      /  V       2N  +  l"C 
whence 

dV N d^ 

V  2  +  (4  —  n)  N  •    C 

For  the  same  range  and  a  different  V  and  C  to  find  the 
change  in  ^: 


See  Appendix  II. 
See  Appendix  III. 
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d_sin2^  2N  +  1      2+(4  — n)N    dV      2N  +  1       N      dC 

sin2^^  N  +  1  2N  +  1       'V       "N  +  1  •2N  +  1  C 

whence 

d  sin  2  ^    _  __  2+(4  — n)N     d^  _     N        d  C 

sin  2^  N  +  1         •    V         N+1  *  ^C 

that  is 

^1        •    o              2+(4  — n)N      -,,      ./        N        ,,      ^ 
d  log-  sm  2  ^  =  — j^— ^^ .  d  log  V  —  ^  ^  ^  .  d  log  C 

These  are  illustrations  only;  other  applications  will  suggest 
themselves. 

THE    SOLUTION    BY    QUADRATURES    IN    THE    CASE    F  (v)    =  A,  v" 

The  equation  of  the  hodograph  becomes 

d  (v  cos  0)        V  F  (v)         A„  v"-'i 


d^ 

gC 

gC 

Placing 

V  cos  d  sec 

<p   =   M, 

we  have 

du 

cos^^^  - 

It^" 

'^  cos"^i  <P  sec 

ni-l^ 

or 

du 

K 
gC 

,  cas"  (p  .  sec"' 

'ddd 

whence,  by  integration  between  V,  (f,  and  u,  d, 

^^  +(^^n=—r~^rP--^~-    +  (f).  -  constant  =  (0„ 


n  A„  u"  cos"  ip  "     n  A„  V"  cos"  (p 

and  the  substitution  of  the  value  of  u  in  terms  of  0  and  ^  enables 

the  problem  to  be  solved  by  quadratures. 

Thus  the  expressions 

-, ,  u  cos  (f  .  .       n 

dt  =  — ^  d  tan  d 

g 

dx  =  u  cos  (p  dt 

dy  =  u  cos  (p  tan  0  dt 

ds  =  u  cos  ip  sec  0  dt 

become  on  representing  ^-^  by  /^"  (which  causes  the  equation 

of  the  hodograph  to  assume  the  form 

(0„-WJ 
d  tan<^ 


\ 

^u  cos  ^ 

dt  == 

k 
g 

dx  = 

g- 

dy  = 

k^ 
g 

ds  = 

g 

dtang 

tan  0  d  tan  9 

(W.— (^).)''" 
sec  6  d  tan  S 
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By  the  use  of  these  formulas  and  by  taking  small  incre- 
ments of  0  so  that  a  mean  value  of  0  will  apply  in  the  incre- 
mental expression,  successive  increments  may  be  computed  and 
added.  The  summation  of  Ay  between  (f  and  zero  gives  the 
maximum  ordinate  and  —  <^  is  that  value  of  0  in  the  descending 
branch  for  which  the  total  summation  of  Ay  is  zero.  It  is 
generally  found  by  interpolation  between  the  last  two  values 
of  d  used. 

The  Quadratic  Law. 

This  process  is  best  illustrated  by  the  use  of  the  quadratic 
law  of  resistance,  as  that  has  practical  application  in  mortar 
fire,  and  in  curved  fire. 

The  equation  of  the  hodograph  becomes 


U^  COS^  (f  V^  COS^  (p 

and  the  expressions  for  dt,  dx,  dy  and  ds  become 

,        k^       d  tan  0 T   C  d  tan  0 

g   '  V(d-iO)  ~        X'^gA^  *  VW—  (0) 
k'      d  tan  6/  C       tan  6>  d  tan  ^ 


dx  = 
dy  = 
ds  = 


g  •   W-W     "       2A.,  •      W  —  (0) 

k^      tan  0  d  tan  0  C       tan  d  d  tan 


g   •       (0  —  (^)  2  A3  •       (0  —  (0) 

k^     sec  0  d  tan  0 C_     sec  <9  djtan^ 

g  '  "  (0—W     ~      2A3  • "  w  — W 
which  may  be  put  into  the  forms 

dt  dtan^ 


^2  g  A3  . 


VC  vld-{d) 


dt' 


oA      dx  ^        dtan^  _   ., 

^A3    .    ~^-  (c)  —  {0)  ^^ 

^  .       dy tan  6  d  tan  6    _    ,  , 

^  A3  .   ^  -  (^^  _  ^^^^      -  ay 

^   .      ^ sec  0  d  tan  6    _      , 

^  A3 .  ^  -  ^^^  _  ^^^      -  as 

It  is  thus  clear  that  dt',  dx',  dy',  ds',  are  dependent  only  on 
i  and  0  and  that  £  is  known  when  -y~^  and  (f  are  known,  since, 
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'  V 

It  is,  then,  possible  for  a  given  ~j^  and  (p  and  d  to  calculate 

t  X 

and  tabulate  t',  x',  y'  and  s'  from  which  —=^  ,  7^,  etc.,  may  be 
calculated  when  A  and  g  are  known. 

Having  calculated  t',  x' ,  s'  etc. ,  for  the  point  of  fall  we  may 

produce  ratios  eliminating  C,  g  and  A2  from  all  except  the 

V 
ratio  -=:.  .     Thus,  for  the  point  of  fall, 

V  C 


rp. 

T'^ 

C 

2gA, 

X 

c 

X' 

2  A2 

s 

S' 

c 

2A2 

and  these  being  divided  by  -7^  give  the  values  of 

T      X       S^ 

The  maximum  ordinate  is  found  by  summing  dy'  from  <p  to 
0  so  that 


c     2A2 

and  this  may  be  divided  by  ^  ,  or  ^^  and  -^^  or  ^  tabulated. 

The  velocity  of  fall  is  found  thus: 
The  equation  of  the  hodograph  gives 

U^co  COS"  <p 


V^  cos^  <P 


(^)  -  {<p) 


Hence 


That  is 


(^)  —  (<p)  _  v"«o  cos^^secV 

1(0  — (^) 
cos  (p  secco  ^  ()4^ 


V       """  ^  \  (0  +  H 

It  is  also  true  that  ds'  is  directly  integrable,  thus 
.  ,  _  _  sec^^d^  _  _      d(^) 
"^^  (0-W  (0  -  (0) 
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Hence  * 

id  -  id)     '^^^  (0  -  (0 

and  this  gives  a  means  of  readily  performing  the  calculation  by 
quadratures,  as  follows: 

1.  Beginning  with  any  desired  -=^  and  (f,  calculate  (0. 

V  c 

2.  Calculate  As'  for  1°  intervals  from  ^  to  0°  and  then 
from  0°  to  beyond  —  ^. 

3.  Designating  the  successive  values  of  0  used  as  0^,  6^  etc. 

n        _i       n 

take  a  mean  value  6'  of  6  as  ^^ — o"^^^  ^^^  ^^^  ^^^'  ^^  ^^^  P^^' 

and  calculate  as  follows: 

Ay'  =  sin  ^'As 
The  summation  from  f  to  0°  of  Ay'  gives  yo',  and  the  sum- 
mation should  continue  until  y'  =  0  in  the  descending  branch, 
at  which  point  is  the  point  of  fall  and  6  =  —  co. 

4.  Ax'  =  cos  6'  As 
Sum  from  ^^^  to  —  <^,  for  X'. 

5.  The  values  of  v  at  each  point  are  found  from  the  equa- 
tion of  the  hodograph  thus  

V0  =  V  cos  <p  sec  6  .   y)        W 
and  a  mean  v  over  the  arc  is  taken  as  the  half  sum  so  that 


^    =  2 


6.  At'  = 

Sum  between  ^  and  —  co. 


Vp  +  Vp+i 

i 
As' 


In  the  Table  of  High  Angle  Fire  Ratios  given  in  Part  I,  the 
quantities  given  are 

V       lOX      lOT  V.      yo       .1       h 

TT '  ^T  '  ^^' '"'  T'    X  ^'^^  ^""^  X 

^  is  found  by  multiplying  ~-  by  —   as  determined  by  using, 
f^  -^  Yo 

in  place  of  (p,  a  mean  of  ^  and  co,   in  the  expression  already 
deduced  for  a  mean  height. 


APPENDIX  I. 

The  Altitude  Factor. 

Consider  two  points,  one  of  them  being  h  feet  above  the 
other.  Let  d,  p,  and  t  represent  the  density,  pressure  and 
temperature  at  the  height  h  and  d^,  p^,  to,  those  at  the  height 
zero,  or  at  the  lower  point. 

Then,  for  equilibrium 

dp  =  -  0^  g  dh  (1) 

and  for  a  gas  p  =  Hag  d  T  (2) 

in  which  T  is  the  absolute  temperature  Fahr.  corresponding  to 
height  h,  or  T  =  461 +  t,    H  the  height  of  the  homogeneous 


atmosphere  at  0°  Fahr.,  and  a  =  ^iy     Hence,  since  H  = 

=  24600, 

we  have  Ha-  53i. 

From  equation  (2) 

dp      do       dT 

p        0^    "^    T 

(3) 

and  from  equation  (1) 

dp  -  —  <?  g  dh 

(4) 

Hence,  from  (2)  and  (4). 

dp              dh 

p            HaT 

and,  from  (3), 

do^           dT         dh 

d  T        HaT 

Next,  assuming  a  drop  of  temperature  of  b°  per  foot  of 
height  (b  =3^0  as  found  by  observations  in  this  country  and 
abroad) ,  we  have 

dT  =  —  b  dh 
and  by  integration  between  T,  and  To,  y  and  0, 

T  -  To  =  -  bh 
or  T  =  To  -  bh 

Hence 

^=  (    1         .V,d(To-bh) 
d       VHab      '■''^     To  — bh 
and  integrating  between  h  =  0  and  h  =  h. 

whence 

d       ^-      V       Tn  / 
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Using,  To  =  520°  Fahr.,  corresponding  to  to  =  59°  Fahr. 

H  a  =  53  J 
b 
this  becomes 


30^ 


^^      \       1560 

.0064. 


(: 


156000/ 

]^         \  —  4.625 

1000  / 

The  value  of  h  for  direct  fire  is  |  yo  and  for  high  angle  fire 
is  a  function  of  the  elevation  and  range  as  given  in  Part  II, 
under  the  discussion  of  the  trajectory  for  which  C  =  oo,  h 
as  used  in  both  these  cases  being  the  mean  height  of  the 
trajectory. 

For  a  mean  height  7500  feet  the  value  of  f^  may  be  repre- 
sented by  /^  =  1  _  AAAAoai^  J  which  form  holds  for  high  angle 

fire  with  sufficient  accuracy. 


APPENDIX  11. 

The   Coefficient  of  Form   as  Determined   by  the  Shape  of  the 
Head  of  the  Projectile. 

Experiments  conducted  by  the  military  authorities  in  this 
country  and  in  France  have  furnished  sufficient  data  to  clearly 
indicate  that  the  coefficient  of  form  is  proportional  to  the  mean 
sine  of  the  angle  made  by  the  tangent  to  the  ogive  with  the 
axis  of  the  projectile.  This  mean  value  is  taken  over  the 
ogival  surface. 

Let  ABC,*  in  the  figure  represent  half  of  the  longitudinal 
section  of  the  ogival  surface.  D  is  its  center.  D  A  =  n  calibers 
is  its  radius.  A  0  =-  i  caliber.  0  Y  and  0  X  are  the  coordinate 
axes.  Let  B  be  any  point  on  the  circle  ABC,  which  curve  by 
revolution  around  the  axis  0  X  generates  the  ogival  surface. 
Let  the  angle  a  =  G  B  F  represent  the  inclination  of  the  differ- 
ential surface  d  S  to  the  axis  0  X. 

The  equation  of  the  ogival  section  is,  then 
(y  +  n  —  i)'  +  x'  =  n' 
and  hence 

dy  2x , 

=  —  ^ -~^ ::  =  —  tan  a 

dx  2y  +  2n^ —  1 

Taking  the  general   case  of  a   surface  of  revolution  the 

generating  curve  being  not  necessarily  a  circle  we  note  that  sin  a 

obtains  over  an  elementary  surface  27ryds  =  27rysec«dx  =  dS. 

The  mean  value  of  sin  a  over  this  surface  is 

1    sin  a  d  S        1  2  TT  y  tan  a  dx         |27rydy 

sm  oc^,  ==  J_o =  J^ =  ^ 

S  S  S 

For  a  complete  head  the  y  limits  from  bourrelet  to  point 
are  i  and  0. 

Hence,  for  a  complete  head 

TT  V"  J  o  ^ 

But  c  =  1  when  the  head  is  a  complete  ogival  struck  with 
a  radius  2;  and  when  c  is  desired  for  any  other  head  we  divide 
sin  a  by  sin  a.^. 
*  See  illustration  on  next  page. 
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Thus 

_  sin  a  _  ^ 
sin  a^       S 
That  is  for  a  continuous  generating  curve  concave  toward 
the  X  axis  at  all  points,  the  coefficient  of  form  is  inversely  pro- 
proportional  to  the  surface  of  the  head. 

PARTICULAR  FORMS  OF  HEAD. 
1.     An  Ogive  Struck  with  a  Radius  n  Calibers, 
In  this  case 


/^ Y      2  n 

y  =  i/  n'       x' 

.  — 1 

2 

He:                 ^^X 

1  n'  —  x' 

dS=  27ryds 

Hence  between  x  =  0  and  x  =  Xi 

s=/%„(.---. 

dx 

V  n'  — x' 

f            2n-l 

;in-i    ^' 

-     ^^  V  -1  2       *  ""^         n 

Now  for  a  complete  head  x^  is  found  by  placing  y  =  0,  when 


Xi  =  t  1    4  n  —  1 
Consequently,  for  a  complete  ogival   head,    n-caliber  radius, 

S„  =  TT  n  (i  /4n-l  -  (2  n  -  1)  sin  -  ^  ^^"~M 
^  2  n       / 

=  TT  n  (  1/  4n  — 1  —  (2  n  —  1)  A) 

in  which  A  is  the  angle  (in  radians)  whose  sine  is ^ 

^  n 

For  n  =  2,  we  have 


2  77  f  1    7    —  3  sin 


^^~) 


g 
The  values  of     (^"  ,  or  c,  for  heads  with  radii  from  2  to  7 

are  given  below : 


n 

c 

2 

1.00 

3 

0.82 

4 

0.71 

n 

c 

5 

0.64 

6 

0.58 

7 

0.54 

These  agree  with  the  most  recent  experiments  at  Indian 
Head  and  at  Sandy  Hook. 
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2,     Ogival  Heads  with  Snubbed  Points, 

Where  the  point  of  a  projectile  is  cut  off  by  a  plane  perpen- 
dicular to  its  axis,  it  is  said  to  have  a  snubbed  point. 

Let  m  be  the  number  of  inches  cut  off  from  the  length  of 

the  projectile  and  represent  this  length  in  calibers;  that  is  -r  , 

by  m'.     Then  the  mean  sine  from  the  base  of  the  head  to  the 
point  Xi,  yi,  i.  e. ,  to  the  flat  nose,  is 

sina= ^it^ZZi!) 

TT  n  (2  Xi  —  (2  n  —  1)  sin  - 1  — ^- 
\  n 

and  for  the  flat  face 

sin  a  =  1 
The  first  of  these  obtains  over  an  aunulus  of  the  cross  sec- 
tion of  the  projectile  whose  area  is  tt  (|  —  y^^),  and  the  second 
over  the  area  of  the  flat  face;  which  area  is  tt  y^^ 

Multiplying  each  by  the  area  over  which  it  obtains  and 
dividing  by  the  sum  of  the  areas  which  is  the  area  of  cross- 
section,  we  find 

sin  .  =  -^-^       -Mi-y^ ^  _^   r^yl 


^(2x,-(2n^l)sin-i-^) 

=  ^  (1  -  4  y,^)^ 

n  (2x1  — (2n-l)sin-i  ^-^) 

In  this  equation, 

Xi  =  i  i/4n— 1  —  m' 

/^^ 2n-l 

1  —  1    n  —  Xi"  "2 

Hence  for  a  snubbed  point 

sin  a       o  00^    • 
c  =  -. =  3.824  sm  a 

sm  ^2 

Representing        sin  ~  ^    ^  by  A 

n  (2  n  —  1)  A  by  A' 
n  (2  Xi  —  (2  n  —  1)  A)  by  B  =  2  n  Xi 
i  (l-4y,^)^byD 
we  have 

c  =  3.824  (-§-  +  4  y,^) 


+  4y,^ 
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Ex.  A  6  inch  projectile  whose  head  is  struck  with  a  radius  of  7 
cahbers  has  its  point  snubbed  (a)  2Y';  (b)  IV;  (c)  f'.  Find 
c  in  each  case. 


Here 


m 
m' 

Xi 

2  n  Xi 
A 
A' 

Yi 
D 
B 

B 

sin  a 
c 


(a) 
2.25 
0.375 
2.22308 
31.12312 

.32318 
29.40938 

.1376 

.21356 
1.71374 

.12405 

.07574 
.19979 
.764 


(b) 
1.50 
0.250 
2.34808 
32.87312 
.34207 
31.12837 
.0946 
.23242 
1.74475 

.13321 

.03580 
.16901 
.646 


(c) 
0.75 
0.125 
2.47308 
34.62312 

.36108 
32.85828 

.0486 

.24530 
1.76484 

.13899 

.00945 
.14844 
.568 


The  general  formula 


c  = 


can  be  applied  to  any  complete  head  of  parabolic  or  other  form 
provided  its  curvature  is  continuous  and  concave  toward  the 
axis  of  the  projectile.  Snubbing  can  be  dealt  with  in  a  manner 
analogous  to  that  used  in  the  case  of  ogival  projectiles. 

The  formula  given  in  Part  I,  viz. : 

2      I  4  n  —  1 

^=^T^ 7— 

assumes  that  c  is  proportional  to  the  value  of  a  at  the  point. 
This  may  be  shown  as  follows: 

2x 


^an  -   -    2  y  . 

f2n  — 

1 

and  y  = 

0  at  the  point, 

and  X 

= 

Hence  at  the  point 

tan  a  ^  y    4 

n  — 1 
n  — 1 

- 

whence 

sin  a  = 

tan  oi 
sec  OL 

= 

sec  a  = 

V 

n'- 

1 

2n 

-x^ 

r  V  4  n 
sec  OL  = 

(2 

n  — 

1) 

V  4n  — 

1 

2n 
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Hence 


sm  oi^  =  — L 


and 

^  _   sin  a:    _  _4^     /  4 n  —  1    _  A    /4n  — 1 

sina2  2n    \        7  n  \        7 

This  formula  gives  sensibly  the  same  results  as  the  formula 

^       S 
but  is  not  applicable  except  to  ogival  heads  which  are  complete. 

Conical  Snubbed  Points. 
These  are  such  that  the  point  is  not  cut  off  flat  but  a  short 
cone  used  instead.     The  inclination  of  the  elements  of  this  cone 
will  lie  between  90°  and  the  point  inclination  of  the  complete 
head. 

As  already  shown,  sin  a  for  the  complete  head  is  given  by 
the  equation 

sina  =  l^4n-f 


2n 
and  sin  90°  =  1. 

Taking  the  inclination  of  the  elements  of  the  short  cone  we 
may  find  c  by  interpolation  between  the  value  of  c  found  for 
the  complete  head  and  that  for  a  snubbed  point. 


APPENDIX  III. 

Experimental  Determinations  of  F  (v) . 

The  value  of  F(v)  has  been  determined  with  considerable 
accuracy  in  recent  years. 

The  values  used  in  our  present  ballistic  tables  are  from  the 
Krupp  firings  at  Meppen  and  the  formulas  through  which  they  are 
applied  are  due  primarily  to  Mayevski  and  finally  to  Zabudski 
who  extended  the  formulas  to  values  of  v  as  high  as  3600  f  .s. 

^rom  current  firings  it  is,  and  has  been  for  some  years,  evi- 
dent to  the  writer  that  these  laws  are  only  approximately  correct; 
and  upon  a  careful  examination  of  the  firings  conducted  by  the 
Gavre  Commission  the  respects  in  which  they  are  inaccurate  are 
clearly  shown. 

At  the  same  time,  certain  interesting  relations  appear  from 
the  firings  of  the  Gavre  Commission.  It  appears  that  Newton's 
Law  of  the  Square  is,  in  the  mean,  true;  since  if  we  take  the  veloci- 
ty 330  metres  per  second  (1083  f.s.)  and  its  experimental  F(v), 
as  determined  by  the  Gavre  Commission,  it  admits  of  the  deter- 
mination of  a  parabolic  relation  y  =  Kv'  giving  a  parabola  asymp- 
totic to  the  actual  curve  F(v)  =y. 

F(v) 
This  is  made  plain  by  plotting  the  curve        2    ="  y'  and  the 

straight  y'  =  K.  The  latter  is  clearly  the  asymptote  of  the  former. 

^  -7  The  law  of  air  resistance  seems  to  be  intimately  connected 

with  the  velocity  of  sound.     As  given  by  the  firings  of  the  Gavre 

Commission  it  almost  obeys  the  law. 

F  (v)  ^       ,      b  (v  —  vj 

v^         ^  "^  c^  +  (V  — v,)^ 

in  which  v^  is  the  velocity  of  sound  in  air. 

Inasmuch  as  the  Krupp  experiments  do  not  agree  with  the 

Gavre  results  altogether  and  the  latter  are  not  capable  of  complete 

expression  by  the  equation  connected  with  the  velocity  of  sound, 

it  is  considered  advisable  to  reproduce  the  Gavre  laws  otherwise. 

F(v) 
Representing  by  K  (v)  the  ratio       2   ,  it  is  noted  from  the 

experiments  that  K  (v  — 1083)  +  K  (v  + 1083)  =  2  K  (1083) .  That 
is,  the  curve  is  central  with  respect  to  the  velocity  of  sound  1083, 
which  is  the  value  indicated  by  the  firings.     The  general  mean 
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velocity  of  sound  deduced  by  Rowland  is  33 1.78  metres  per  second 
at  0°  C,  or  1089  f.s.,  in  pure  dry  air. 

The  minimum  value  of  K(v)  occurs  when  v  =  330 — 155  m.s. 
or  v=175  m.s.  =  574  f.s.,  and  the  maximum  value  for  v  =  330 
+  155  m.s.  =  485  m.s.  =  1591  f.s.;  the  mean  of  which  is  1082.5 

f.s.     The    corresponding    values    of    K(v)    are    -jTrr  and    -.L 

o  0^7  o  on 

giving  a  mean  of  ^^r.     The  value  ^(yr  is  given  for  1083  f.s.  by 

the  experiments.     These  facts  suffice  to  determine  the  curve,  if 
it  is  of  the  form  outlined  above,  since  in  that  case 
a  =  .0000887 
b  =  .0000933c 
c  =  508 

This  determination  gives  a  curve  astonishingly  like  the  experi- 
mental curve  and  leads  to  the  supposition  that  some  curve  of  that 
form  must  represent  the  law.  Maxima  and  minima  are  hard  to 
locate  on  an  experimental  curve  and  the  probable  error  of  the 
determinations  is  sufficiently  great  to  render  it  possible  to  find  a 
curve  of  the  above  form  giving  correct  results. 

The  probable  cause  of  the  errors  of  determination  which  un- 
doubtedly exist,  is  the  failure  to  record  the  exact  condition  of  the 
atmosphere  with  respect  to  the  ratio  of  the  specific  heats,  which 
exercises  a  marked  infiuence  on  the  velocity  of  sound.  For  this 
reason  the  result  is  more  or  less  composite  in  the  vicinity  of  the 
velocity  of  sound,  which  in  ordinary  air  at  mean  temperature  and 
humidity  is  higher  than  1089  f.s. 

In  our  present  tables,  given  in  Artillery  Circular  M,  the  laws 
used  were  those  of  Zabudski,  based  on  the  Krupp  firings.  They 
are  of  the  form 

F  (V)  =  A,  v-" 
over  certain  velocity-ranges,  and  rendered  compatible  by  the 
adjustment  of  A  when  n  is  changed.  In  integrating,  to  find  the 
S,  A,  I,  and  T  functions,  constants  of  integration  are  of  course 
necessary  to  ensure  compatibility  in  the  determination  of  these 
functions.  (See  Appendix,  Artillery  Circular  N.)  These  velocity 
ranges  are  sometimes  quite  extended  and  the  mean  value  A  as- 
sumed as  a  coefficient  of  v"  departs  materially  from  many  of  its 
variable  values;  ^side  from  this,  the  determinations  are  not  the 
most  recent  ones. 

It  is  thought  best  to  give  the  value  of  K(v)  as  found  from  the 
firings  of  the  Gavre  Commission  and  quoted  in  full  by  Charbonnier 
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(Traite  de  Balistique  Exterieure,  1905) .     The  use  of  these  values 
in  computing  the  ballistic  functions  will  appear. 

In  Table  I  the  values  of  v  are  taken  at  such  intervals  that 
ordinary  or  linear  interpolation  only  is  required  between  tabular 
values. 


TABLE 

I. 

Values  of  10^  K  (v)  = 

F(v) 
V- 

X  10' 

V 

10'  X  K  (v) 

V 

10'  X  K  (v) 

f.  s. 

f.  s. 

0 

0.5385 

1214 

1.1784 

263 

.4689 

1247 

1.2246 

492 

.4202 

1279 

1.2581 

558 

.4136 

1312 

1.2862 

590 

.4133 

1345 

1.3071 

623 

.4147 

1378 

1.3213 

656 

.4172 

1411 

1.3357 

689 

.4206 

1443 

1.3452 

722 

.4276 

1476 

1.3524 

787 

.4484 

1509 

1.3558 

820 

.4658 

1542 

1.3580 

853 

.4866 

1575 

1.3594 

886 

.5149 

1607 

1.3561 

919 

.5481 

1673 

1.3527 

951 

.5902 

1903 

1.3038 

984 

.6465 

2297 

1.1995 

1017 

.7168 

2526 

1.1505 

1050 

.7967 

2920 

1.0879 

1083 

.8863 

3150 

1.0637 

1115 

.9761 

3510 

1.0445 

1148 

1.0561 

3609 

1.0411 

1181 

1.1262 

3937 

1.0325 

1214 

1.1784 

4000 

^   1.0308 

Table  I  gives  the  rneans  of  interpolating  so  as  to  find  K(v)  for 
each  f  .s.  from  0  to  4000  f  .s.  Multiplying  in  each  case  by  v^  we 
obtain  F(v).  We  may  thus  readily  compute  by  simple  quad- 
ratures the  S,  A,  I,  and  T  functions,  adding  increments  from  one 
velocity  to  the  next  at  intervals  of  1  f  .s.,  or  more. 

To  do  this  beginnings  with  4000  f.  s.  the  value  of  K(v)  for 
which  is  10.308,  we  have 


d  S  (u)  = 
d  I  (u)  = 
d  T  (u)  -= 
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u  d  u 


F  (u) 
2gdu 

u  F  (u) 

F  (u) 

dA  (u)  =  I  (u)  dS  (u). 
The  S  and  T  functions  may  be  made  equal  to  zero  at  4000 
f .  s. ;  also  the  A  function.     The  initial  value  of  the  I  function 
may  be  determined  as  follows: 

Represent  K  (v)  by  K,  for  simplicity. 

Placing  F  (u)  =  K  u^  =  -^^^ —  u%  we  have 

dI(u)=-2-?A^  =  -^^ 
u  F  u  K     u"^ 


Hence 


2  ^  1  S" 

I  (u)  =  -^  X  „— 2~  =  =r^  +  constant 
'  ^         K         2\f        Ku^ 


Ku 


+  k 


Now 


Hence 


d  A  (u)  =  I  u  d  S  u 

i^ju ki^     du 

K'u'        K  •    u 


A  (u)  =  2^^  ~  K"  ^^^'  ^  +  k« 


If  now  ki  be  made  zero  and  k^  =  —   ov^uomV    ^^  follows 
that  ki  being  zero 

^  ^^^^^^  ^     K  (4000)^  ^  "FT4^000r" 
The  initial  values  having  been  determined  the  solution  is 
one  of  ordinary  computation  by  quadratures. 


UF  820  .H3  pt.2 

Hamilton,  Alston,  1871-1937 • 

Ballistics  • • • 


MAR 


NE  CORPS  U  LIBRARY 


3000044518 


Library  of  the  Marine  Corps 

2040  Broadway  Street 
Quantico,  VA  22134-5107 


iOBi  i^tBm 


-a^ 


